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ON SEVERAL ASPECTS AND APPLICATIONS 
OF THE MULTIGRID METHOD FOR SOLVING PARTIAL 
DIFFERENTIAL EQUATIONS 

By Nathan Dinar* 

College of William & Mary 
Williamsburg, Virginia 23185 

Research Sponsored by NASA Langley Research Center 
Contract NAS1-14972 -3 


SUMMARY 


Several aspects of multigrid methods are briefly described 
in this report. The main subjects include the development of 
very efficient multigrid algorithms for systems of elliptic 
equations (Cauchy-Riemann, Stokes, Navier-Stokes) , as well as 
the development of control and prediction 'tools (based on local 
mode Fourier analysis) , used to analyze, check and improve these 
algorithms. 

Preliminary research on multigrid algorithms for time 
dependent parabolic equations is also described. This report 
deals also with improvements in existing multigrid processes 
and algorithms for elliptic equations. 

Some partial and typical results are given. More complete 
and detailed information will be presented in the author's 
Ph.D. Thesis, to appear at the Weizmann Institute of Science, 
Rehovot, Israel, 


I . INTRODUCTION 

This report deals with several aspects concerning multigrid 
methods for fast solution of partial differential equations. It 
covers the research on this subject for the period August 1977 - 
August 1978, when the author was spending his sabbatical at the 
NASA, Langley Research Center. This research is part of a Ph.D. 


*Research Assistant, Department of Mathematics 
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Thesis to appear shortly at the Weizmann Institute of Science, 
Rehovot, Israel, including more detailed results and conclusions. 

The research on multigrid methods began in the early 1*970' s 
by the initiative of Professor A. Brandt. He is today supervising 
several projects on these subjects, including the present research. 

The multigrid method can be applied to a wide range of 
problems and, therefore, it interests many people. It is known 
to be one of the most powerful and advanced methods used today. 

The multigrid method uses the fact that the numerical 
discrete equations we usually want to solve are not independent. 
They are derived from a continuous problem whose solution we 
want to approximate. In the process of the solution of the 
discrete equations it is convenient to keep in mind the differ- 
ential origin of the problem. The use of discrete operators 
on several levels of meshes, interacting strongly in the process 
of the solution, allows us to solve the problems on the finest 
grid very efficiently, with a minimum number of arithmetical 
operations (0(N)), where N is the number of equations in the 
the finest grid. 

The multigrid method consists generally of several 
processes, performed in a given order, and defining an iter- 
ative cycle. These processes include generally: 

Relaxation: (Usually of Gauss-Seidel type) Used only to 


smooth the errors, i.e. to reduce these high 
frequencies in the error that are not described 
in coarser grids. 
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Transfer of Residuals to a Coarse Grid : This allows us to 

define a problem on a coarse grid, that is 
similar to the original one. The solution on 
the coarse grid is used to improve the approxi- 
mate solution given on the finer grid. 
Interpolation : Used to define a new approximation on a 

finer grid, given an approximation solution on 
a coarse grid. 

Moreover, we can intermix the principles of the multi- 
grid method with the grid adaptivity principles, which mean 
adaption of the order of discrete approximations and mesh 
size, using total smoothness of the solution. This keeps 
N as small as possible. 

These basic ideas and others, as well as treatment of 

theoretical and practical aspects have been extensively 

12 3 

covered in the papers of Brandt ' ' ' , 

This report deals with a- wide spectrum of problems 
involved in the development of multigrid methods and multi- 
grid software and the principal parts are the following: 

a) Development of multigrid methods and software 
for elliptic systems of equations. (Including 
Cauchy-Riemann equation, Stokes equation, and 
Wavier- Stokes equation. 

b) Development of control and prediction tools for 
various steps and processes involved in the 
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multigrid method. (Based on local mode Fourier 
analysis) . 

c) Development of multigrid methods for time 
dependent parabolic equations. (Heat equation 
on rectangular domains) . 

d) Improvement and changes in existing multigrid 
processes and algorithms (mainly on Poisson 
equation) and treatment of new ideas in relaxation 
methods . 

II. MULTIGRID METHODS FOR SYSTEMS OF EQUATIONS 
The basic ideas of the multigrid method are not 
restricted, of course, to a unique equation, and from a 
theoretical point of view no special problem could be 
expected in implementing multigrid ideas to a system of 
equations. However, special and detailed algorithms for 
this problem did not exist, and it was extremely important 
to get sharp and practical proofs of the efficiency of 
multigrid methods for systems of equations. One of the 
questions we did not know the answer in advance, was for 
instance, what is the appropriate method of relaxation for 
a system? 

In order to answer this and many other important 
questions, we developed multigrid algorithms for three 
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typical problems, starting with a simple one, and each 

subsequent problem involving more and new complications, 

relative to the former. Moreover, since it is usually 

wise and necessary, in this kind of research, to isolate 

different questions which may arise, we developed the 

algorithm for simple geometrical domains (rectangles) . 

The past experience showed that more complicated geometries 

did not affect the efficiency of multigrid methods (see 

4 

for instance Shiftan ) . For complicated geometries there 
are, of course, more programing problems. 


2 . 1 Cauchy-Riemann Equations 
The equations are: 

V v y = 

V V y * f 3 (x ' y) 

for 0*x<l; 0<y-M. 

And the boundary conditions : 

U (0,y) = U Q (y); v(x,0) = V Q (x) 
u(l,y) = u 1 (y); V(x,l) = v^x) 


Figure 1. 
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(3) 
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The data must fulfill an additional condition, due to 

the continuity equation (2), in order to assure, a solution 

to the system. This condition is 
i // 

](U 1 -U 0 )dy+ J (V-j-VgJdx = J j f 3 (x,y)dxdy (4) 

o 0 - «■ 

The discretization of the problem was accomplished 

on a staggered grid. This method has several advantages 

in problems of this type. For instance, it allows us to 

define easily second order conservative finite differences. 

The grid is described in figure 2a. 

Figure 2a. | j j 



The finite differences approximations for equation - (1) 
are defined at grid intersections (X), for instance: - 

t[ U 3- U l- (V 2-V 1 >j = fi (5) 

And the finite differences for equation C2) are defined in 
the middle of the cells (.), i.e.: 



(6) 


As in the continuous case, the data must fill an additional 
discrete condition equivalent to (4) in the finest grid, i.e.: 
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k >U-V+ L 2 < V li- V Oi> = Z < 4a > 

J i' J 

Several relaxation schemes appropriate to this system 

were considered and all these schemes have a good smoothing 
factor*, around -^=.5. The most convenient scheme was chosen 
because of its simplicity and its small number of arithmetical 
operations. This method belongs to a new approach of relax- 
ations, developed by Professor A. Brandt, and is called 
"Distributed Relaxations". 

This approach is characterized by the fact that when 
passing through a point (or cell) in the grid, for which 
the difference equations are defined, we change the value 
of more than one unknown, in order to make zero residual 
on this point. 

In the Cauchy-Riemann equations separate relaxations 
are performed to each one of the equations, according to 
the schemes discribed in figures 3 and 4. 


*Efficiency of relaxation m a multigrid process is measured 
by the "smoothing factor," which is defined as the asymptotic 
ratio between the fast components of the errors after a relax- 
ation sweep, and the same errors before the relaxation. A more 

2 

exact definition can be found in Brandt . More considerations 
on this subject can be found in Chapter 3 of this report. 
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Figure 3 ; Relaxation of U x -V^ = 


V, to 


u 2 +r 




U 1 -f 


Figure 4 : Relaxation of U x +V y = ^3 

! V 2 + f' 


U-. i f. 
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The values U+J and V+ c in figures 3 and 4 represents the 

•*> 

new values after the sweep of this point. <? is always chosen 
to make the residual zero at the point. The relaxation is a 
Gauss-Seidel type relaxation and it passes through all the points 
.(or cells) of the grid in a usual natural order. The important 
property of this method relies on the fact that when relaxing 
each of the equations, the residual of the other equation 
remains unchanged . 

Mode Fourier Analysis shows that the smoothing factor for 
this method is Jj -.5, the same as in Poisson equation, which 
is, of course, equivalent to the Cauchy-Riemann system. The 
multigrid algorithm shows convergence factor -*j X .56* 
which is about the same value as in the Poisson case. The 
number of operations in the relaxation of the Cauchy-Riemann 


*See Chapter 3 in this report 
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system is, of course, higher than in Poisson equations (about 
two times more). On the other hand, we can say that the 
information we get from the solution of the Cauchy-Riemann 
system (functions U and V) is also two times the information 
we get from the solution of the Poisson equation. 

The main reason for our treatment of Cauchy-Riemann 
system was for learning purposes. It is probably the simplest 
elliptic system possible and it allows us to treat compli- 
cated problems more easily. However, this system is inter - 
estxng in itself, as we can see in the works of Lomax and 
Ghil^, These papers describe fast algorithms for the 
solution of the Cauchy-Riemann system. We think that the 
algorithm just described in this report is preferable, 
because it does not depend, in principle, on the simplicity 
of the domain, and because it is at least as fast as Ghil 
algorithm, and even faster asymptotically (for finer and 
finer discretizations) . The implementation of the algorithms 

for nonlinear equation is also much easier in the multigrid 
3 

method . 

An example of a computer output is given in the Appendix, 
Output No . 1 : 


2.2 Stokes Equations 
The Stokes equation are given by 


Au-p = f. 

(7) 

x 1 

A v-P y = f 2 

(8) 


U +V 
x y 


( 9 ) 
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The unknowns are U, V, and P. U and V are given on 

boundaries. Like in Cauchy-Riemann systems. The data must 
satisfy the condition (4) . The discretization is also in the 
same staggered grids already described. P and f 3 are defined 
in the centers of the cells, f^ and f^ are defined in the 
same points as U and V, respectively (see figure 2b) . 



The difference equations for equation (7) are given by 

Au(x,y)--)p(x+k ,y)-P(x- b- ,y ) 1 = f , (x,y) (10) 

k, i^L 2- J, j ± 

for (x,y) where U is defined. A is the usual five 

7- I - J 

points discrete Laplace operator: j 1 ~‘ i ‘ 

lx L i J 

Similarly, for equation (8) 

AjVfxry)- ijp(x,y+ )-P(x,y- b )J = f 2 (x,y) (11) 

for (x,y) where V is defined. 

For the continuity equation (9) , the difference 
equations were already defined by (6) . 

Equations (10) and (11) hold for interior points. Near 
the boundary the difference equations are different, but they 
are simple and will not be described here. 
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The relaxations are based on the same principles already 
described. The relaxation of equations (10) and (11) , the 
Momentum equations, are performed separately, by the usual 
Gauss-Seidel method corresponding to the Poisson operator. 

The relaxation scheme for the continuity equation (6), is 
described in figure 4. 



$ is chosen so that equation (6) is fully satisfied on the 
current cell, and the changes of the values of the P are 
chosen in a way that preserves the residuals of equations 
(10) and (11) when relaxing equation (6).* 

The theoretical smoothing factor in the relaxation 
method just described can be shown to be *j =.5, like in 
Poisson and Cauchy-Riemann equations. 

In the numerical experiment we get a multigrid conver- 
gence factor** -*j%.65» This is a very good value, and assures 

*This is exact only on unbounded domains. On bounded domains 
this property cannot be strictly maintained near the boundary. 

**This is the reduction factor of the residuals per unit work 
which is equivalent to one relaxation sweep on the finest grid. 
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a very fast method for solving the Stokes equations. However, 
it was a little worse than our expectations, based on the 
previous experience of Poisson and Cauchy-Riemann equations. 

In order to check this point, and to reduce the possi- 
bility of errors in the computer program, we performed a 
complete multigrid Mode Fourier Ana-lysis (with the aid of the 
techniques to be described later) . The results of this analysis 
showed complete agreement between numerical and theoretical 
results, leading to the conclusion that the value of ^ £ .65 
is intrinsic to this system, which is characterized by strong 
interaction between the equations. The Mode Fourier Analysis 
helped to rule out possible programing errors or bad influences 
of the boundary on interior regions. 

An example of a computer output is given in the Appendix, 
Output No. 2: 


2.3 Navier-Stokes Equations 
Two-dimensional Navier-Stokes equations are given by 




(12) 

O-aK+v^) 

]v- P, - ■O.'j) 

(13) 

4- 

- 3 -- 

(14) 


and the same boundary condition as in Stokes equation. The 
parameter R is called the Reynolds number (R=0 corresponds 
to the Stokes equation) . 
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These equations are more interesting, from the practical 

and physical point of view, than the former equations. The 

numerical treatment is more difficult and complicated, due 

to several reasons, like the nonlinearity and the fact that 

boundary-layers may appear for large Reynolds numbers. 

These facts demand a careful choice of the difference 

equations, in order to keep the elliptic ity of the difference 

operator. This can be accomplished in several ways, based 

usually on one-sided first differences, instead of central 

first differences, that may cause instability, unless we make 

a very drastic and unpractical restriction on the mesh size h 

1 

(h must be of the order h=0 (^-) ) . In this work, we define 

the following difference approximations to the first derivatives 

of U and V, for example for ^ — we define 


* TC C 


(15) 


for O^a^l (a x = 

a is defined by 
x 


,5 corresponds to central differences) 


and ^ = k-RU (x,y) 




X 


|j)4l 

1 


1 3( 


- — 


14- 


Vx r 


and 




(16) 


are defined in 

a similar way. It can be shown that with this approximations 


we get finite differences with the desired properties. The 
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equations, defined on the same staggered grid already 
described, are given by 



t i> 

- P(x- 



(18) 

V^r 

t L ?(V 


r-j] 5 








(19) 

jc/cv^ i0 ) 

-UK-b 

,» v 

tV'u,y 


[y, 3 b 

(20) 










i 


-s ft?) 


(21) 












y^kRv ■ h = 

kRV 

(22) 



• ” 


, . 1 

P w 




^ i 
' x . 

| + (/- 


1 

(23) 


■ 5 = 

<f(j>+ fr-j) + t fty + fr-'i) 


(24) 


{These are the equations for interior points. The operator 
is a little different near boundaries) . 

The relaxation method for (18) and (19) are done as in the 
Stokes equation, by freezing U and V in the calculation of 
(f(i) and ft?)- 


The relaxation scheme for equation (20) is described in 
figure 5. 

original pagejs 

Of POOR QUALM 
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in the current cell is zero. The changes in P would keep 
the residuals on (18) and (19) unchanged if we froze the 
coefficients. 

Mode Fourier Analysis shows that the smoothing 
factor (calculated for equations with frozen coefficients) 
will generally depend on the direction of the relaxation of the 
Momentum equations. For instance, if the solutions U and V 
satisfy U?0; V^O, and the direction of the relaxation is as 
usual, (say by columns and increasing y, and increasing x) 
then Ajz ,5. If U or V are negative in some part of the 
domain, and R is big enough, the smoothing factor will be 
higher, and, therefore, the multigrid method less efficient. 
These difficulties may be easily overcome by making relaxations 
sweeps in different directions, improving in this way the 
smoothing factor, or better, by the technique of distributed 
relaxations already mentioned. 

The first preliminary version of the algorithm does not 
include the features just pointed out, and, therefore, it works 
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with the full multigrid efficiency only for cases where U and 
V do not change sign in the domain, or for general cases with 
R not too big (up to R=100) , 

These first experiments are very important in showing 
how the algorithms work in principle independent of secondary 
problems, like the direction of the relaxations, mainly 
technical. 

The numerical results shows a good multigrid convergence 
factor, bounded for all R by ^ ^ .7-. 75 . This provide* 

a very fast algorithm for solving Navier-Stokes , and 
almost the same efficiency as in the linear Stokes case. 

i 

The full multigrid cycle (explained later) was also 
implemented. Preliminary experiments shows that only 
8-9 works units* are needed to solve the equations 
to the level of truncations errors. (Execution time 
on CDC 6400 for a grid of 64x64, is about 15 seconds^ 
and the program is not optimized, and includes all kind 
of calculations for debugging purposes) . An example 
of a computer output is given in the Appendix, Output 
No. 3. 

III. CONTROL AND PREDICTION TECHNIQUES 
IN MULTIGRID METHODS 

One of the remarkable advantages of the multigrid 
method is the possibility of predicting in advance its 
efficiency, given the relaxation method, interpolation, 


*Equivalent to relaxations sweeps on the finest grid. 
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etc. This can be done by means of local Fourier Mode 

Analysis in an infinite space, i.e. neglecting boundaries 

influence. This can be justified because of the fact that 

Fourier Analysis is very accurate in descibing the fast 

components behavior, while is less accurate for slower 

components. However, most of the calculating work in the 

multigrid processes is invested in the reduction of the 

fast components of the errors in the finest grid, and this 

is achieved by means of relaxation. More detailed justi- 
i / 

3 

fications can be found in Brandt , who shows a very good 
agreement between theoretical results based on Fourier Mode 
Analysis and between numerical data obtained by multigrid 
algorithms. 

We can distinguish two levels of analysis that are 
appropriate for multigrid methods. In the simplest one 
we calculate the smoothing factor , which depends only 
on the relaxation process. The more complete and complicated 

analysis takes into account all the multigrid processes and 

0 ^ 
estimates theoretically the value of ^ {definitions of Sf 

o 

and are given in the footnotes in Chapter 2) . In many 

cases, it is enough to perform the simpler analysis. In 

.0 - 

these cases, the formula % ~ , where f is the 

hvrs*. 

ratio of the number of points between coarse approximate 

and fine grid, (usually J-i), and d is the dimension of 

o 

the problem provides a fair approximation to ^ 
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The calculation of the smoothing rate involves the 
search for extremum of functions in bounded domains, and the 
calculation- of eigenvalues (defined usually by a general- 
ized eigenvalue problem Ax=ABx) of complex matrices of 
size <2 >.<2 , where ‘j is the number of equations in the 

system. The estimation of the convergence factor 

involves the same kinds of calculations, but the size 

cf ci 

of the matrices are much bigger ( 2 x 1 ^ where a. 
is the dimension of the problem).* 

It is clear that the kind of calculations just 
described, cannot be performed in a closed form, 
with the exception of extremely simple cases, like the 
smoothing rate for Poisson equations and one dimension 
simple problems. 

Because of this and other important reasons, an 
algorithm that performs these calculations was constructed. 

A general FORTRAN subroutine was written for this purpose, 
and the main inputs are the following: 
d - The dimension of the problem. 

<3. - The number of differential equations = the 

number of unknowns. 

r - The number of relaxations in one multigrid 
iterative cycle. 


*In the relaxation process, the Fourier components are 
independent. But when we take into account the transfer 
from grid h to grid 2h we get interdependence between 

groups of 2 U components. 
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A,B - Complex matrices that are derived from the 
particular difference equations and methods 
of relaxations. The size of these matrices 
is ^ and they .depend on the Fourier 

component defined by 

A ^ _ A function of the Fourier component & , 
u 

describing the initial distribution of 
amplitudes in the errors. (For a random 

initial distribution of errors, A^ =1 
for all $ ) . 

L - Number of points defining a mesh in the & domain. 

f — The ratio between the number of points on grid 

and the number of points on grid 
The output contains mainly: 

a) The smoothing factor /y defined by 


where 








dtb [A (&)- 


and J t$l - Av-'O, x j y; { 

r 4 ' ><i 

b) The weighted smoothing factor given 

the number of relaxations r and the weight 


e * 


( 25 ) 


( 26 ) 


function A 
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c) The location (mode) where the smoothing 
factor '‘j is found. 

d) A distribution map of (only for the 

two-dimensional cases', d=2') . 

The same subroutine, with some changes, can be used 
for the complete Fourier Multigrid Analysis. In this case, 
the user has considerably more work in defining the complex 
matrices A and B that include, m this case, the transfer 
of residuals from grid h to grid 2h, and the interpolation 
from grid 2h 'to grid h. 

These programs are very important in the research of 
multigrid methods. It can be used for several purposes, 
for example: 

1) Checking new relaxation methods or new multigrid 
processes, before the algorithm is translated into 
the computer. 

2) Comparison of several algorithms for the purpose 
of optimization, i.e. looking for the faster and 
stable ones, 

3) Debugging of multigrid computer programs. 

We used the subroutine in several cases, including 
Cauchy-Riemann and Stokes equation. We also used it for the 
Poisson equation. In all these cases, we did the complete 
Fourier multigrid analysis, checking and comparing various 
methods of relaxation, residual weighting and interpolation. 
We got full agreement between theoretical and numerical 
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results (in all cases considered) , and this fact increases, 
of course, the reliability and prediction power of this 
theory . 

An example of these facts are summarized on table 1, which 
include theoretical and numerical results for Stokes equation: 
Table 1 


Comparison Between Theoretical and Numerical Results 
_ fpr Stokes Equation ' ’ 


*- 

(The number of 
relaxations on 
the finest grid 
in one iterative 
multigrrd cycle) 

— 

r* 

Parameter related 
to the degree of 
accuracy in solving 
coarse grid correct- 
ion equations 

'Y 

Multigrid Convergence Factor 

Theoretical multi- 
grid complete 
Fourier analysis 

— 

■Numerical 

results 

1 

.4 

.661 

. 638 

1 

.5 

.648 

.634 

2 

.3 

.680 

. 695 

2 

.4 

,710 

.725 

3 

.1 

! 

4 

.714 

.722 


The subroutine was also applied for checking a special approach 
to distributed relaxations, that may perhaps be applied to 
Navier-Stokes equations. This will be shortly described in 
Chapter 5. An example of a computer output is given in the 
Appendix , Output 4 . 


IV. MULTIGRID METHODS FOR TIME DEPENDENT 
PARABOLIC EQUATIONS 

A possible quite obvious way to use multigrid procedures 
for initial value problems is considered in the work of 
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3 

Brandt . By this procedure we use a multigrid algorithm 
for solving the implicit equations usually defined at each 
time step. If, for instance, we want to solve the Heat 
equation in two space dimensions then for each time step, 
an elliptic problem similar to the discrete Poisson equation 
is defined (and, in this case, the Gauss-Seidel relaxation 
has even better smoothing properties than in the Poisson 
case) , The typical amount of work needed in advancing each 
time step by multigrid procedures will be, accordingly, 
equivalent to 5-6 relaxations (see ^ ) . If a solution for 
M time steps is required,’ then the total amount of work 
will be about 6M' relaxations . 

The question which naturally arise is whether we can 
use multigrid principles and ideas to get more efficient 
methods for these problems. 

The answer appears to be positive, and Professor A. 

Brandt pointed on some approaches which eventually can develop 
in efficient multigrid algorithms for these problems. 

A basic idea is that marching in time can be done for 
most of the time steps on coarser grids. 

The appropriate equations on a coarse grid are carefully 
corrected in a way that assures the correct representation 
of the information from the fine grid, so that even when 
marching on coarser grid, the accuracy of the finer grid is 
kept. To update the information from the finer grid we need 
some sporadic and infrequent time steps on the finest grid. 
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. that constitute, of course, most of the computational work. 

i 

The base of this approach relies on the fact that fast 
Fourier components of the solution (represented on finer 
grids) converge to steady-state after a very short time. 

Changes in the solution after this are due mainly to slower 
components, that also change slower in time. This allows 
us to march on coarser grids with large time steps, after 
the influence of fast components have disappeared. 

In order to check these and other ideas, we developed 
two different algorithms for the Heat equation, in a 
rectangular domain with Dirichlet Boundary Conditions. The 
first algorithm uses the Crank-Nicholson implicit scheme, 
and the second one uses the simplest explicit scheme. It 
must be pointed out, that the stability restriction in the 

in- 
explicit scheme that makes this method very unpractical, 

does not appear in this case, because most of the time we 

march on very coarse grids, where large time steps are allowed. 

a. The Implicit Scheme 

The equations to be solved on the finest grid 
(with mesh size h^ 1 and time step k°) are: 

- ?(*,%) (27) 

and the unknowns are u^ (x,y,t+k°) for each (x,y) defined 
on the grid. The index 0 in u^,h^,k^ represents marching on 

* * 

j5r 4 i 

k* H 
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the finest grid. Similarly u m ,h ra ,k m will represent marching 
on a grid which is m levels coarser. The equations on coarser 
grids are given by 


~ jA /W ^v^t)4 U)3;+T ^( Xl . V| S) 

, "" c 

where usually k- -1 H and l represents the appropriate 
’ rt 

correction. X" has a very important significance, and it 
represents the spacial truncation error of grid m, relative 


to grid 0 . T. 


is given by 


rv 

/Vl \ 

T k ° - T (A -4, )a ‘"k jji 

^ Si 

V 

where 5 - (C /; ^ } satisfying A\ ^ c > «• 

The transfer from a given coarse grid m to a finer grid 
m -1 is done by ■ 


m -1 m -1 , _m-i m m -1 . 

u = "Last + X m - "Last 


X 

where r ~ means interpolation (cubic) from coarse 
to fine grid and "Last" represent the last marching in 
time on the m -1 grid. 

b. The Explicit Scheme 

The equations on the finer grid are: 

r ■ 

~ l *£(* ,* 5 .t + k") - -U f (X.'f) 

L ” . K J 


o? *°° 
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On coarser grids the following equations are defined: 


- A r fcty + T^(X, j, S) 


( 31 ) 


In this algorithm we keep a constant ratio A over 
all grids given by A~ ♦ The value of A - .2 

was chosen as very appropriate from the point of 
view of the fast Fourier components. 'L represents 
here the complete relative truncation error (in 
space and time). u xs defined as follows: 

For m=l 




L 




and 


S, < t 


Then for general m 


T 


i ^ 

A 


M 

-I 

j 5 -’ 


i-i 


r ■ 

Jv 1 




(32a) 


and s is defined as in (29) . The transfer from coarse 
to> finer grids is performed lake in (29a) . 

It must be pointed out that this research, which deals 
with new approaches in implementing raultigrid ideas, has 
a very preliminary character and the principal aim was 
to check, potentially, the various promising possibilities. 
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In this context, some numerical experiments were performed, 
in which we tried some criteria in order to define the 
times where we switch from coarse to finer grids and vice- 
versa, In addition, we compared the accuracy we got in these 
algorithms with the same schemes when we marched in time on 
the finest grid only. The numerical results we got from 
these experiments, in both algorithms, are good and encouraging. 
However, the explicit scheme appears to be preferable in some 
aspects. The numerical results in these experiments, and some 

theoretical considerations arising from the interpretation 

» 

of these results, lead to the conclusion that applying these 
multigrid ideas with full efficiency, (i.e. solving the problem 
in an amount of work comparable with the work invested in 
solving elliptic problems) , means using adaptive techniques 
where we can change and adapt the order of the discrete 
approximation . 

But even without adaptive techniques we can use the 
present algorithms to solve the Heat equation very efficiently. 
The amount of work needed will depend on the smoothness of the 
solution. For instance, for a problem with smooth initial 
conditions, we can save 98% of the work, in relation to the 
same numerical scheme defined on the finest grid only (of 
size 64x64) . 



V. IMPROVEMENTS AND CHANGES IN EXISTING 
MULTI GRID ALGORITHMS 
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In this chapter, we deal with existing multigrid proce- 
dures, and we check various approaches in order to improve 
their efficiency. The first part deals with the practical 
solution of Poisson's equation in a small number of numer- 
ical operations. Although in this part we concentrate on 
Poisson’s equation, for reason of convenience when checking 
and comparing numerical results, it is quite clear that the 
techniques developed here can be applied to other elliptic 


problems as well. 

The second part of this chapter deals with a special 
approach to distributed relaxations for general five points 


difference operators of the form 


a 


d 

-s c 


s = a+b+c+d. 




b 


.1 


5,1 The Poisson Equation 

The first numerical multigird algorithm and program 
was written for the Poisson equation. For this equation 
there exists today more information related to multigrid 
methods than for any other problem. 

In Brandt' 1 ', the algorithm called there Cycle "C" 
for solving the Poisson equation is described with detail. 
In Cycle "C" the multigrid iterative cycle starts at the 
finest grid, where an initial approximation to the solution 
is defined. This approximation is then improved by solving 
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correction equations on coarser grids. The switching from 
fine to coarse grid and vice versa is controlled by some 
internal, criteria'. 

Cycle "C" is very useful in learning and understanding 

multigrid performance and several theoretical aspects, 

like the asymptotic convergence factor, etc. It is 

not generally the most efficient algorithm for solving 

real problems, for which we do not know in advance a’ 

» 

good approximation on the finest grid. In addition, it 
is usually difficult to know in advance how many iterative 
cycles are needed to get the desired accuracy. One can 
perform several iterative cycles and then get very close 
to the solution of the difference equations but this will 
generally be wasteful, because in a real problem, we do 
not need more accuracy than the accuracy defined by the 
differential truncation error . 

Because of these and other considerations a fixed 

' * , 3 

algorithm for Poisson equation is described m Brandt . 

In this algorithm we start with an approximation on the 
coarsest grid, and after we perform one iterative cycle 
on each level,' the solution on each grid is used as a first 
approximation on the next finer grid, by means of a .cubic 
interpolation. This algorithm is based on the knowledge 
and past experience in solving Poisson equation by multi- 
grid techniques, and provides a solution to the problem 
up to an accuracy comparable with the truncation errors 
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accuracy, and this solution is found in a minimal number 

of arithmetic operations (we do not even need to calculate 

the residual norms, used generally' for. internal criteria) . 

In the present work a similar but more general algorithm 

is constructed. This algorithm can serve for learning 

purposes as well as for solving practical problems. It can 

also be easily applied to other elliptic operators. The 

algorithm was implemented in the FAS (Full Approximation 

Storage) mode. (For a detailed description of FAS see 
2 

Brandt ) . In this mode we can look at the coarse grid 
correction equations in a somehow dif ferent, * but equivalent 


way, as follows: 


, Ik 

L *v v 5 V " * 


2J. 

tH is a new and better approximation than to the 

k 

exact solution /^ n of the difference equation 




(34) 


IL < 

given by (33) describes the truncation error of the 

grid 2h, relative to grid h. It is well 'known that the exact 

solution u of the differential equation Lu~f satisfy the 

following 
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L 2h U = £+ ~2h 


(35) 


and X 2 ^ is the differentia-1 truncation error of the 2h 
grid. The comparison between (35) and (33), that arises 
naturally during the PAS Multigrid Cycle without additonal 
investment of computational work, shows the remarkable 
differential aspect of the multigrid method. This is in 
sharp contrast with other methods, where the algebraic 
aspect is usually the most important and central one. 

These properities open new and interesting possibilities 

u ^ 

in implementing the FAS algorihm: L j in (33) satisfies 

the following 




L k 


(36) 


and the truncation errors can be represented sby ~ } 

(provided additional smoothing conditions) . Then we have 




T ik- l l k + 0(U) 


(37) 


Then we can change the FAS algorithm by defining the 
following correction equation on a coarse grid (instead 
of 33). 



(/ 


lk 


f + i 

3 


T 


'i / 


(38) 
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A series of numerical experiments were performed to 
check these interesting points, including printouts of ij~~E ■'!/_ 

and I I ^ II, ■ (we chose for this purpose problems where 

L/ ?o 

the exact differential solution u is known) . As a result 
of' these experiments, we can point out the following 
conclusions: 

a) The assumption that after only one multigrid cycle 
(5-6 work units) we get the needed accuracy in the 
solution (truncation error accuracy) was confirmed. 
In addition, truncation error behaves clearly like 
0 (h 2 ) . 

I 

b) We can get the same desired accuracy without even 
performing a complete multigrid cycle, i.e. we can 
stop the process before we interpolate back to finer 

A 

grids. The solution is then defined on coarse grid 
points only, but the accuracy is the finest grid 
accuracy . This feature can be very helpful in 
the development of multigrid methods that use small 
amount of storage, much smaller than the number of 
unknowns in the finest grid. 

c) Using thd T-extrapolation (38 instead of 33) we can 
get a much better approximation to the solution of 
the differential equation on the finest grid, much 
better in fact than the approximation achieved by 
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the exact solution of the difference equations 
in this grid, where the accuracy is bounded by 
truncation errors. In order to get this better 
approximation, we do not have to invest any addi- 
tional work I Our numerical results do not show, 
however, that the approximation is 0(h 4 ), in 
spite of (37). This is because of the fact that 
at least 5-th order interpolation is needed for 
this purpose. 

Nevertheless, the improvement we can get by this 
small change in the FAJ3 algorithm is impressive. 
The improvement will generally depend on the 
smoothness of the solution. The more smoothness 
the better approximations. 

In table 2 we can, compare the results we get in the 
different variations we tried. The problem is Au =* f; 


the solution is known, given by u = sin (3x+3y); (0^x<3j 0sy£2) 

Table 2 


Size of the 
Finest Grid 

|| u K _ u\l u ^ 7 

Experiment 

1 

Experiment 

2 

Experiment 

3 

Experiment 

4 

16 x 24 

1. 4xl0** 2 (23) 

1.8xl0" 2 (5.9) 

5 . 3x10*” 3 (5.9) 

2.3xl0~ 3 (9.1) 

32 x 48 

3 . 4xl0~ 3 (23) 

4.2xl0~ 3 (5.6) 

4.7xl0” 4 (5.5) 

2 . 3xl0 -4 (8.9) 

64 x 96 

8. 6xl0“ 4 (23) 

l.OxlO” 3 (5.4) 

5 , 2xl0” 5 (5.4) 

1.8xl0“ 5 (8.9) 
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Key to table 2 : 
Experiment 1 - 

Experiment 2 - 


A/ , 


^ is the exact solution of the difference 
equations (up to randoff errors) . 

v ^ 

0 is the approximation after only one cycle 
in all grids. 


Experiment 3 - 0 A is the approximation after only one cycle 

in all grids, but using (38) instead of (33). 
Experiment 4 - Like experiment 3, but the number of relaxation 

in each level was doubled. 

The numbers in parenthesis represent the number of work 
units which are equivalent to the total number of relaxations 
on the finest grid (printed in each line of the table) . 


5.2 Distributed Relaxation 

One of the important application of the method of distributed 
relaxations is in finding relaxation schemes with good smoothing 
properties. This is generally hard for very asymmetric operators. 
In these cases, the usual Gauss-Seidel relaxation may be very good 
or very bad, depending on the direction in which the relaxation 
is performed. The problem is even more complicated in case of 
non-linear difference operators, because in this case we may need 
one direction for some parts of the domain, and an opposite 
direction for other parts. In these cases we can perform, 
as we explained in a previous chapter, a series of alter- 
nating direction relaxation sweeps, but this cannot generally 
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be an ideal solution. If we think, for example, of a three- 

3 

dimensional problem, we need there at least 2 =8 relaxations 
on each level for a multigrid' iterative cycle, in order to do 
relaxations in all possible directions. This may be too much, 
because we know that in most multigrid procedures, we use 
only 2-3 relaxations for achieving optimal multigrid efficiency. 
So, alternating direction relaxations may in some cases reduce 
considerably the efficiency and, in addition, it will cause 
more programming difficulties. 

Because of this and other reasons it is convenient to 
find relaxations, with good smoothing properties, that are 
to some extent independent on the direction of the relaxation, 
so we can relax always in the same direction, no matter how 
the operator behaves. 

As a first example, we consider the one-dimensional 
operator represented by 1 - ij / 1 1 ^ j r C < ^ <' oo 

which can be derived, for example, from the differential 
equation V £«./- o (fry taking one-sided first differences). 

V nj 

The smoothing factor can be easily calculated in this case, 
and for extreme values of r) it behaves in the following way: 


/ V*v 


1 




1 



I 0 


Relaxation from left to right. 


Relaxation from right to left. 


( 39 ) 


for *1 * o , ^ s . $ 

In the present work, we check a special principle for 
the construction of distributed relaxations. The relaxation 
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is always' done from left to right according to the scheme 
represented in figure 6. 

Figure 6 . 

Relaxation 
■ 




y ,+ *i 

^ »i 






t- 1 




£ and in figure 6 represent the values of the 

approximation before relaxing on point i, and the new values 
after the relaxation on this point are represented by / 

. J is chosen so that the residual at point i is 
zero. The changes performed at points i and i-l make the 
residuals at the points i-1 and i-2 nonzero, although they 
were previously zero when we relaxed the-se points. So we 
can choose the parameter < so that the deterioration of 

these residuals is minimal, say in the L„ norm. 

" z 

If we chose in this way it can be shown that the 
relaxation scheme is always stable and satisfies ^ t .7 

for all o . For . The optimum parameter 

depends, of course, on but satisfies ,4£ ^ .5 

and, in fact, we can chose a constant suboptimal £* , say 

= .45, without significant lose of efficiency. 

It must be pointed out that this scheme is not 

I 

necessarily the best one for all values of ^ . In fact, 

7 

for large rt , there exists better stable schemes (Brandt ) . 
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Using the same principles just stated, we calculated 
distributed relaxation schemes for general five points 
operators in" the form N 

d 

a -s c 

b 

✓ 

where a,b,c,d>0 and a+b+c+d=s. 

Operators of this type are typical, for instance, 
in Navier-Stokes problem. The relaxation scheme is 
described in figure 7. 

Figure 7. f 


> * 




-ifr 


Mrt* 


Direction of the 
relaxation. 


« rfcc' 


oL and jb are the solutions of an appropriate optmization 
problem. As in the one-dimensional case, M + T r t\+oLc ,pi+ft>T 
represent the new values after the relaxation of the central 
point of figure 7 . 

The results (calculated by the algorithm described in 
Chapter 3) show that the smoothing factor is bounded far 
below 1, even in the cases of great asymmetry of the operator, 
except in the cases where the problem degenerate practically 
to a one-dimensional problem (like the case aic c~ 0). 

In the case of Poisson's equation (a=b=c=d=l) 315, which 
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is significantly better than in the normal Gauss-Seidel 
relaxation (^ = .5). For this Poisson operation we get £*=/>-—, * 
For other nonsymmetrical operators, the comparison between 
this distributed relaxation and the usual Gauss-Seidel is 
even more remarkable, as we can see on table 3. 

Table 3 


Operator : Smoothing Factor 


a 

b 

c 

i - ~ 

d 

! 

Gauss-Seidel 

Relaxations 

Distributed 

Relaxations 

.5 

1.5 

1.5 

.5 • 

i 

.631 

.313 

10.5 

.5 

10.5 

j 

.5 ! 

.913 

.449 

.5 

100.5 

100.5 

.5 ; 

.990 

.481 

.5 

1000.5 

1000.5 

i * 5 

.999 

i 

.484 

.5 

.5 

1000.5 

j i 

j 1000.5 jj 

1.000 j 

.752 


Like in the one-dimensional case, it may be possible to find 
more efficient scheme for the extremely asymmetric cases. But 
the present methods assures a good smoothing factor in all the 
cases, bounded far below 1, and consequently, an efficient multi- 
grid algorithm. Moreover, if the problem is nonlinear, we do 
not have to know in advance the behavior of the solution, in 
order to define an appropriate relaxation scheme. In addition, 

*This distributed relaxation scheme for Poisson’s equation 
was implemented in a multigrid computer program, and we get 
perfect agreement with the theoretical smoothing factor 
Jj * .315. 
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as in the one-dimensional case, it is not necessary to invest 
extra work in the exact determination of the parameters and 
They can be estimated in a suboptimal way, without significant 
loss of efficiency. 

It must be pointed out that in the case of the Poisson 
equation the Gauss-Seidel relaxation can be still considered 
a little more efficient, if we take into account the number of 
arithmetic operations needed for each relaxation method. On 
the other hand, for all other operators, even rf they are 
small perturbations of the Laplace operator, the proposed 
distributed relaxations are much more efficient than the 
Gauss-Seidel relaxation. 
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.1 M'S -0 7 


,7b«lE-nR 

.im-07 


j 1 0 h F • 0 7 

3 R , 90 9 

,7l2E*ua 



£-l)b _- 

R3qp.nR 


^U7F-l)7 


..cn 

•EMf) CYCLE WO, ’<> 


--END CYCLF^NO, |0 _ . _ 


•ewd cycle wo. n 


's 


■£wp CYCLE wo. 12. 



j 3 

35^125. . 

.52«;p-(ifl" 

.29«E*U« 

.UJOE-UG 

. HbOE-O* 


.35.929 

. 1 63E-0* . 

.1 t tE-IIH 

. I6HE-H3 . 

.21 7E-OK 

! 5 

35.911 

.10NE-IH 

.7 s 0 e - 0 *» 

• l?"E-OA ... 

. IMF -9 9 

.3 

„ 3b.^»u« 

,2b5E-'J6 

• 17NE-0?. 

.2R0E-06. 

.3«9F.-»h 


1 1 

.675£-fib 

.96KF.||fl 

.7b?E-'lR ... . 

.XIHF-Ort 

5. 

_ . 3*.?bl 

. 1 *) 1 E • 0 7 

. Q 37c*"y 

. U2E-07 ... 

. 2 0 3 b - n 7 

6 

_ ,37.261 


.utdE-i'7 

. . .hme-i)7__ 

. 9 ' * E " 0 7 

h 

__ J«.26l 

.'i7«E-tJ7 

.2"7E*n7 

,i(SiJE-(i7 

.7«5F-b7 


■end cycle **o. 13 


i — 


Cif CUE 

NO. 

2 

CYCLE 

EFF. 

.691 

. ACC. EFF. 

• o«! . .... 


_ _ , . a ^ 

CYCLE 

NO. 

3 

CYCLE 

EFE. 

,520. 

ACC.FFF. 

.563 



CYCLE 

NO. 

'1 

CYCLE 

fcFF. 

.£65 

ACC.FFF. 

,b«3 



CYCLE 

Nn. 

5 

C YFI E 

FFF. 

.5 aft 

ACC. EFF. 




CYCt.F 

Ml. 

6 

CYCLE 

fff. 


ACC.FFF. 

.S91 



CYCLE 

Nil, 

7 

cycle 

EFF, 

.6*0 

a C C , f E F , 

.597 



CYCLE 

NO, 

a 

CYCLE 

EFF, 

.636 

ACC. EFF. 

.-'02 . . .. 



CYCLE 

Nfl, 

9 

CYCI E 

FFF. 

.635 

*CC. EFF. 

.606 



CYCLE 

«n. 

1« 

CYC> E 

EFF. 

. 6'l2 _ 

... ACC.FFF. 

, b 1 0 


... 

C YT| f 

Ml, 

1 1 

CYCLE 

fff. 

.65? 

ACC. EFF. 

.FlO 



CYfl.F. 

o.n, 

12 

CYCLE 

EFh. 

.656 

ACC.FFF, 

, t> 1 9 



CYCLE 

on. 

1 3 

CYCLE 

EFF . 

, 665 

ACC.FFF. 

.621 

- 

• 


“ 

' ■ 

CovviAoeucB Factor f 

„ co 

■— 

- “• 




FOR GYlCl* CVCLc 

. . ' . FA CTnIV . 
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- 
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. ~ 

— 

. . - 

— 
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--- 

s. 

— 

— - ‘ - -- 

■ -- 

— 


— 
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- 

- 
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OUTPUT No 3 


H^rtMHM-HKMrfMhMWMHrtWHW-l-'Mrtt'khl^HWHkMWHrlMWMriHWWWHHMwMWHwMWMWHhHWMViMWHWMWMNMWHWWWMhMWHWMWMWHSMNM-rtWM^HWMWMiNMHrtHh^^"^”^"^^ _ __ 

^ U ^K« rt «K!-MIV«hM|.^y»MWHHK*MM 1 ulMHMWhWk'Ut>HMWMWrt»)>tWhWHlNMW>'WM l <VkHkMWHWt<VtHWMWMWUWMWHW»<WMWMWMWMWMWMWHWMHMkMV(K'<HWMt»HHMV(M l HMWM»lV.kMHH«MH ^ 


EQUATIONS 


NAVIER STOXeTI 1 1 ) (DEL-HUO/nX-RVO/OY) U«PX«F l 
(2) (DEl-HUIVOV-RVn/OYJV • PYbF 2 -- (3) llX + VY»F3 


BOUNDARY CONDITIONS AMD XMIU+L FUNCTION VALUES -• - 
METHOD - 

2 Y-INTEHVAIS, Mn» 
,300 


- , Ft»F2«F3=0, — -- — 
U*Y*SQRT (H ) *FXY V«X*SttMT(R)*FXY Ps-R**2/2 CX*+2+Y**2) (FXY« 1+X { I »X ) Y ( J »Y 5 ) 
FASrLIKF INJECTION, SFPARATE RF L A X A T I ONS , 1 1 NE A R I N T 


COARSEST GRID IS 
MlU'hER OF GRIDS 
METHOD PARAMETFR5# 
TUL* 
FI). 1 ! ,350 

REYNOLDS NU-MFR 


? X-INTFRVALS, AND 
6 

ET *» ,?fm OELTAa 
,10nE«nu 
,3Rt) ,300 

« q 1 0 E ♦ 0 « 


,500 


LcVEL C = CifXi H G-RiJ) 



ACt. 

kORK 

| 

R«Tio>/iu 

VoiJMJ 

• “ 


LEVEL 

J WElDHTeO 

Six. 1 

Ea X 

sa 3 

_ 

6 

1,000 

, 1 9 5 E ♦ 0 * 

,3705+OR 

, 103E + OS 

,KMF+03 


b 

2 ,noo 

.76RE+07 

.139E+08 

.73AE+U7 

.371E+03 


s 

2,250 

,6535+06 

, 1 «ot+o7 

-- , «7 IE+06 

. ,u7HE+02 


h 

3,250 

,33nE+oo 

, 665E ♦ 0 6 - 

.27BE+06 

- .lluE +03 


b 

R.250 

, 1 il 7 E ♦ n 5 

,3S9E*65 

.616E+0R 

.693E+H? 


b 

5,250 

.U16F+03 

.9URE+03 

.19BE+03 

,5?5E+02 


\ 

* b 

6,250 

^ , 1 60E+03 

' ,256fc + (>3 

. 169F+03 

, H30F + P? 

act k 

5 

6,5.10 

,5i3r.*oi 

,918F +n? 

- ,1365+09 

, 1 <1 ? t + 0 ? 


_. 5 • 

6.750 

- ,?3hF + n3 

. , 8 ONE + 1*3 

- .25RE+03 

, 9R5E + 0 1 


<i 

6,613 

- .31AF+03 

,1 10E+O3 

,7976+03 

,20 IF +01 


(J 

6 , *75 

. . ,2335+03 

, 3>9 2 E + o 3 - 

.27JE+03 

, 1 57E + 0 1 


3 

6,001 

,1535+03 

,725^+02 

,36«t+03 

, J99E+0 0 


3 

6.006 

, 156F+03 

.267E+03 

, 1 RiE+03 

, 1 A2E+00 


2 

6,010 

,5706+0? 

.. , 3 5 7 E + 0 2 - 

, 127E+03 

, U p 1 F »0 1 


? 

J 

2 

? 

6 |ri« 

,71 95. + 0? 

- - , J29E + 03 — 

— ,766E+d? 

- .151E-01 


6 R 1 5 

, 1 27E + 0? 

.. , t 99E + 02 - 

, 170E+02 

, ?9iie»0 2 

6,019 

,5265+03 

.79SE+03 - 

,7rt6F+03 

, 1 06E-P! 


b’,°?\ 

.917F+02 

, 1 56E+03 

. 106E + 03 

,759F»02 

2 

1 

6,927 

6.926 

.77 1E+02 

.937E+01 

. . .126E + 03- 

. 2 0 o E + o 2 

— .9U1E+02 

- .6905 +0-i_ 

.299E-02 
..ouuE-ol • 


fcNO CYCLE NO, 1 


■END CYCLE NO, -2 



2 

6,932 

,1065+03 

,1575+03 - - 

.196E+03 

.235E-0P 

2 

6.936 

,2105+02 - 

,9395+02 -- 

,1665+02 * 

.176E-02 

3 

6.0*51 

, ) 895 + 1)8 

,2905 + 09 

,1635+09 

,1015+00 

3 

6. 

, 1 OftE + oS 

, 1 A*E + P 3 

, 1155+03 

.903E-01 

3 

6,9*2 

,2*>>5+03 

,fl?uE+0? 

.9026+02 

, ? 1 up »0 1 

tt 

T.0OS 

, 81135+00 

, 1 UbE + 05 

.8395+09 

, 992E+00 

u 

7.1U7 

,2065+03 

,9675+03 

,1205+03 

,3*05+00 

a 

7.170 

, 50 o 5*02 . 

. ,0205 + 0? 

, 6 0 6 1 + 0 2 

.2075+00 

5 

7,920 - 

.270F+05 

,5105*05 

,2795+05 

, 68 0 E + 0 ! 

5 

7.670 

,6| 0E+O3 

, 1505+09 

, 198E+03 

.2796+01 

5 

7 ,9?(i 

,6275+02 

,1095+03 - - 

.693F+02 

.155E+01 

5 

8,170 

,3j6fc+C2 

,59V.+n? 

,3526+02 

.101E+01 

b ■ 

0.170 

.5615+05 

.10 35 + 1)6 

, 5 7 3 E + 0 5 

.255E+02 

b 

10,170 

. 1 A?5 + oa 

,952£+09 

,6795+03 

, 1 2 0 E + 0 ? 

6 

11,170 

,6375+02 

, 1 ? 7 5. ♦ o 3 - 

,9*75+02 

.77PF+01 

h 

12.170 

,3965+02 

,69?t+0? 

,2976+02 

,5806+01 

m mm mm 

5 

12.920 

,661) 2 + 0? 

, 1 A65. + 02 

,1796+03 

,1886+01 

5 

12,670 

.339E+U? 

,6*95 + 0? - 

,2725+0? 

, 1 2 9 E + 0 1 

9 

12.73? 

,3175+02 

,2095+02 

,6996+0? 

,3725+00 

a 

12. 70S 

,11 095 + 0 2 

,7065+02 

,9625+02 

,2095+00° 

3 

12, Ml 

,1695+02 

,2015+02 

, 2 6 9 E + o 2 

, 4 9 A F - 0 1 

3 

1 2.A26 

, it 6 1 F * 0 2 

,5805*02 

.738E+0? 

.2355-01 

2 

1 2, *30 

,1075+02 

,16*5+02 

j 1 3«E+02 

,U9«E»02 

3 

12.896 

,2675+03 

.9035+03 

.359E+03 

,1655-01 

3 

12.861 

,5605+02 

. 1205 + 03 

, 9 ? 5 5 + 0 3 

.109E-01 

2 

1 2, *63 

,3005+01 

, 5 3 ?t ♦ o l 

.951E+61 

.228F-02 

3 

12, PHI 

,9335+02 

,1315+03 

, 7 A 1 F + » 2 

,63*5-0? 

3 

12,696 

, 1 3 1 &+ 0 2 

.235E+0? 

,1 3*E+02 

,8 075-02 

3 

12.912 

,5615+01 

, 60 HE + o 1 

, 993f+0] ' 

.3085-0? 

u 

12.975 

.A825+P3 

, 12?5+0« 

, 122E+09 

,7266-01 

<i 

13.037 - - 

,U«Q5+n2 - 

.{(ilE+03 

.2765+02 

, 37*5-01 

u 

13,100 

, 12OK + 02 

,17*5+02 - 

. 1 7 Ht + 0 2 

.2325-01 

3 

13,115 

,9005+01 

,6175+01 

, 53oE+P 1 

,5975-02 

3 , 

13.131 

,SSfi5+ol 

, 7U9E+01 

,799E + 1>1 

.37*6-02 

2 

1 3.135 

, 1 555*0 l 

,3<-«5 + ol 

.1355+01 

, 10*5.0? 

3 

13,150 

,2535+0? 

,?6HE+0? 

,9595+02 

.2325-02 

3 

13.166 

,9605+01 

, 1 97E+02 

.8795+01 

,1 AOF-0? 

2 - 

13.170 

,9b55+on 

,1915+01 - 

, 1 oot+01 

,5976-03 

3 

13,106 

, J«0t+o? 

,?S*++n2 

,26*5+02 

, 1 1 * E -02 

3 

13,201 

.699F + M 

, lU0t+0P 

, 5 9 2 E + 1) 1 

.9715-03 

2 

13,205 

,9'IOF + II u 

, t?71+oi 

. 199F+0I 

,31 55-' l 3 ’ 

3 

1 3.221 

,5135+0! 

, 6rt3E+o 1 

.7*1 E+01 

,6*95-03 

3 

13.236 

,2995+01 

.579E+M 

,2b55+01 

.6356-03 . 

u . 

1 3,299 

, fl u 3 5 ♦ o 2 — - 

. , 1 1 t.E + «J — - 

- , 1 2SE+03 

- , 1 0 6 6-0 1 

a 

13.361 

,6605+01 

, 1 H6E+02 

,5986+01 

,5915-02 

5 

13,611 

, 3965+09 

,9995+09- - 

,9956+09 

,9876+00 

5 

1 3. fi M 

,1 26F+03 

.3325+03 

,2*7E+02 

.295E+00 

5 

19,111 

, 1 0 ? t + o 2 - 

, ! 6 7 K ♦ n 2 

, 123E+02 

.21 15+00 

h 

15,111 

, l 065 +05 

. 1 59E + 05 

.1936+05 

.3556 +0 1 



UD 



S— 

— nrrm 

— ryp^CTua — 

• x JOCTTHT- 

intf+nr 

P27A7T01 

6 

17,111 

,15'iE + n? 

.274F+0? 

.1386*02 

,1808+01 

6 

IMI1 

,8666+01 

- ,t«Hl+02 

• ,4826+01 

, 160F+0I 

5 

ta.361 

, l 89E+02 

- ,4756+01 

- , 8866+02 

.599F+00 

5 

1 fl.bt t 

, 105F+O2 

,23SE+n2 

,62aE+ot 

- .420E+IJ0 

a 

1 8,67a 

, 8 2 U E ♦ o 1 

, a 1 5 E ♦ a l 

- ,i7?E+oa 

. 1 1 2E + 00 

a 

18, 73b 

,997F + 0 1 

, |4«F. + n2’ 

• - .902E + 01 • 

,51 IE-ill 

3 - 

ia,7s? - 

.asae+nj 

, 32 9t +o l 

- -.39AE+01 

.909E-02 

a 

in.aju 

, 190E*03 

, 2308+03 

, 3 t 3 F ♦ 11 3 

• ,?64F-01 

a 

. tB,877 - 

,66»t+0 1 

- ,1558+02 

- -- , 3 5 7 F + 0 1 

- - , 1 4 5 E « 0 1 

a 

18,919 

,?d?F*m 

- , 353E+0 1 

— ,??3F+»1 

,9338-8? 

5 

19, 169 

,9b*F +03 

, i 33 t+oa 

- - , laat +ya 

1 , 1UOE + 00 

5 

19.439 „ 

,23«E+02 

, 622F.+ 02 

-- ,4588+0) 

, 79 | t-0 1 

5 

19,bH9 -- 

, 1 99E+01 

- .3486+01 

- , 1 29E + M 

,5036-0 1 

A 

pn.AflP . 

, 349E + na • 

- .42AE + 04 

,a9bF + oa 

■- *• • 7^6F + 00 

6 

2 1 , 0**7 

,1 19E + 03 

,3138+03 

- - .272E + 02 

.9158+00 

b 

?2,bb9 

, 2 7 6 F + n ] 

,biinE+nl 

- • ■ . I43E + 0I 

.525F+00 

6 

23.989 

*0 1 

,465r + M 

, 8 9 0 F + 0 0 

,0578+00 

5 

23,939 

.417F+01 

, ) 2 9 E + (i l 

.lose +oa 

, 164E + 00 

5 

24,189 

,27RF + (ll 

— ,AaaF+oi 

- -■ , 1 4 | F +0 1 

— , 1 0 9 E + 0 0 - 

U 

?U,252 

,1756+0) 

.1378+01 

,3Mifc+fM 

.259F-0 1 

a 

24,314 

.231 E+Ol 

,4676+01 

• ■ , 19UE+01 

.113E-0I 

3 

?4 , 330 

,84AE + no - 

-- , 6 a 5 E + [i o 

. 120t+ni 

.21 ?E-f'2 

a 

24,393 

,339^*02 

.4K5E+02 

- ,6166+02 

,8188-02 

a 

2a, alb _ - 

.317F+AJ 

. 5 3 3 E + 0 1 

, 3 7 3 F + (i 1 

,38 IE-0? 


a 24 , 5 1 rt .HtiCtOI ,?13t+M ,1778+01 - .24UF-02 


3 

24,513 

, 32»F+no 

,8558+00 

,062E + ||» 

,ba66-03 

a 

?4.596 

. 50*fe+»1 

, 5? JE + 0 1 — - 

,9198+0) 

, I45E-02 

a 

aa.osfl — 

,2)88+01 

, 1746*0 1 — 

,2496+01 

, 1 1 BF -02 

3 

?a,b7a . . 

.282E+O0 

,50 38* + 0 0 

, 3 0 1 F + 0 0 

,3968-03 

a 

24.716 - 

, 3 9 5 F ♦ 0 1 

,« 1 5L + n t - 

.712F+01 

,703F-f>3 

. a — 

. 2a,799 

, 170F + 0 1 - 

,317E +01 

,1848 + 01 - — 

, 706F-03 

3 • 

2 4 , k | a 

,22?F+nn 

,4018+00 - 

,2336+00 

.275F-03 

. £1 

24,877 . 

,32of +0 | 

, 0 a 7 8 + 0 1 

.4678*01 

,5368 -03 

a 

?a,9J0 

. 143F+M 

,2626+01 

.145E+01 • 

,504E-n3 

3 

28,955 

, 1 9 1 E + 0 o 

, 3 2 9 F. + 0 0 - 

,216f +00 

.2058-03 

a 

25,oi« 

,22?E+oi 

,3476+01 

,2688+01 

,8036-03 

a 

25,080 - 

.I09E+1I1 

- , l 99 E + 0 | -- 

,1118+01 — - 

.385E-03 - 

3 

25,090 

, 1 0 1 F ♦ 0 n 

,261 E+OO 

, I 9 8 8 + 00 

.15HF-03 

a 

25,158 

, 1 076 + Ot 

, 23SE + o 1 

, 1 8 5 F + 0 1 - 

, 3 1 3E-03 

a 

2b,??1 

, fl 1 aF + fio 

, 1 87E+01 

,6558+00 

, 30?F-03 

S 

25,471 - 

,2318+03 - 

,3298+03 

,3316+03 

,3996-01 

5 

25.721 

, a on 6 + oi 

,938E + n 1 

.20UE+01 

, 1 9 a F - ( 1 1 

5 

25.971 - - 

,6996+00 

- - , 138E + 01 

,6406 + 00 

, I 2 1 £ - 0 1 

5 

?6,221 

,2HUE + flij 

,5338+00 

,2716+00 

s p?PF.O? 

8 

27,221 

,7018*03 

, ) 07E + O4 - 

, 1 1 OP + 04 

, 1 6 l E +00 

b 

28.221 

, 1 9afc+02 

,5135+02 

.anfiF +01 

,1318+00 

b 

29,221 

.H3JF+00 

- , l 60E + 0 1 - - 

.6H3E+00 

,1078+00 

b 

30,221 

,4858+00 

,9788+00 - 

.3356+00 

.8978-1'! 




ry fxj ^4 \*a c ^ 1/1 *J1 t <> > O 9 * J1 tn r> J1 t/i _/l C fc M fc t; c ^ w Aj ui Aj ^ t^-t 


307*71 ' 

, 66 7 EVAS - 

- ,265E+oo 

".Tme+oT - 

.295E-01 

30,721 

, 5 o 3 E + o o 

, 1 ISE + ol 

- - , 27UE+00 - 

, 1B3E-01 

30.7H3 . 

, 306E+00 

- - , 2?6E + 00- 

• ~ . , 6OOE + 00 

, <1 00E-P2 

30,6<I6 

, 360E+O0 

.625E + 00 - 

- , 0 H3F + 00 

, ifioe-02 

30,861 

- JOflF + OO 

. 165E + 00 

.259F+00 

.332E-03 

30.677 

, .3AIEM10 

.U37E+00 

- .65OF+00 

, 1 7 3 F - 0 3 

30,601 

, 1 J6E + 00 

, 107E + 00 

, 1 U3E.+00 

, 366t-0« 

. 30.60', 

,aoAF+oo 

,6b9k+00 

- - , 0 7 6 £ + 0 0 

,?ii6F-nu 

30.0/16 

,t 17E+00 

,1 AfiE + oO 

- .107E+00 

, 5 0 0 1 - 0 5 

30,690 

,S(iaf + oo 

.659E+00 

- .793E+00 

, 1 OOE-OO 

30.09U 

, 1 ? (1 F + 0 0 

.?()) E + 00 

, 15OE+O0 

, 6 3 0 E - 0 5 

30. 895 

,?26E-»t 

.569E-01 

.771E-02 

, 5 ft l E - 0 h 

30,696 

, 130E+00 

, 1 16E+O0 

.265E+00 

, <1 1 6 F - 0 5 

30. 9(1? 

,3755-01 

,7568-01 

.316F.0) 

.250E-05 

30,9\H 

, 1 035+01 

.2U6E+01 

, 1 6 1 E +0 I 

, 7fl(lE -tl« 

30.930 

,6225+00 

. 1 1 ? E + o 1 

,659t+00 

,5ii6E-ou 

30.936 

,6nlF-ni 

. 1 18E + 0 0 

, 5 0 1 E • 0 1 - 

, 1 5 7 E - 0 0 

30,953 

,5735+00 

, 7 9 1 E + 0 0 

- - ,flU6E+00 

, 3 0 7 F -Oil 

30.969 

, 1 7 l p + 00 

, 3 3 1 E + (i 0 

, 1 5 u E ♦ 0 0 - 

, 179E-00 

3 u . 9 73 

, 107E-01 

, 253E-0 1 

.165E-01 

.U59E-05 

30,928 

, 1 72F+O0 

,20hE+oO 

.205F+00 

, 107E-O9 

3t ,006 

,50i)E»0t 

, ) (iHF +00 

,0655-01 

, 860F-P5 

. 31.(166 

,5 not + oi 

, 6 6 7 E ♦ 0 1 

.762E+01 

.770F-OJ 

31.129 

, 66 1 E + 0(1 

. ISIE + Ot 

— ,9506+00 

,0 17E-03 

31,101 . 

,2125+00 

, <J«eE + 0» 

, 1 60F+00 

,2U6F-03 

. 31,207 

,2106-01 

, 191E-01 - 

, 02 1 E -0 1 

.527E-OU 

31,270 

,mie+o n 

- ,7yuf+00 
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